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Abstract 

The complex Ginzburg-Lzmdau equation with a thermal noise term is studied imder 
conditions when the system is convectively unstable. Under these conditions, the noise 
is selectively and spatially amplihed ^ving rise to a noise-sustained structure. Analytical 
results, applicable to a wide range of physical systems, are derived for the variance, and 
the coefficients and thermal noise term are determined for Taylor-Couette flow with an 
axial through-flow. Comparison is made to recent experiments. 
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Consider the equilibrium state of some spatially extended system and a small spatially 
localized perturbation about this state. If the perturbation grows at a fixed location 
in space, the equilibrium state is absolutely unstable. However, if the perturbation is 
convected with the mean flow such that it grows only in a moving frame of reference, 
eventually damping at any fixed location, the equilibrium state is convectively unstable.^ 
Although little attention has been given to this distinction until fairly recently, it is an 
important distinction since these two types of instability give rise to qualitatively very 
different behavior. In a convectively unstable system external noise is selectively and 
spatially amplified giving rise to spatially growing waves and a noise-sustained structure, 
a concept introduced with studies of the complex Ginzburg-Landau equation. In contrast, 
in an absolutely unstable system structure is sustained by the internal dynaxnics. 

Since noise - whether thermal or otherwise - is an element common to all physical 
systems and since any system with nonzero group velocity will be convectively unstable 
siiificiently close to and above onset of the instability,^’® one would expect noise-sustained 
structure to be very common in nature. For example, in addition to being important 
in classic open-flow fluid systems such as jets, wakes, and channel flow,^“®>’^ the above 
concepts are important in such diverse systems as film flow,®’® binary fluid convection,®’^® 
sidebranching in dendrites, and traffic flow.^ Considering the general nature of these 
concepts, there are undoubtedly systems in other fields to which they would also apply. 

Since the complex Ginzburg-Landau (CGL) equation is a generic equation which de- 
scribes systems near onset of an instability, it has proved to be an ideal system in which 
to explore these concepts. Also, since the equation is rather simple in form, there is some 
hope for deriving analytic results. Considering the fact that, up to this point, no analytic 
results have existed for convectively unstable systems in the presence of spatially extended 
noise, analytic results should prove very useful and would provide further insight into the 
interaction of noise with convectively unstable systems. 
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Recently, noise-sustained structure has been studied experimentally in Taylor-Couette 
flow with an imposed axial through-flow. It was found that under convectively unsta- 
ble conditions, noise-sustained structures of traveling vortices exist. This system is very 
useful for study in that, for sufficiently small Reynolds numbers, noise-sustained structure 
exists in a parameter regime where the flow is laminar and axisymmetric, thus allowing for 
a great deal of experimental control. Further, the CGL equation is valid for this system in 
a parameter regime of experimental interest. 

Since in a convectively unstable system noise is amplified exponentially in space, even 
extremely low levels of noise will be sufficient to produce a noise-sustained structure, 
assuming the system is sufficiently long.^”^ A very interesting question a^ked in ref. 15 is 
whether or not the noise-sustained structure seen in the Taylor-Couette experiments is of 
thermal origin. This is an intriguing and important question since an affirmative answer 
would imply that further efforts at noise reduction - short of decreasing the temperature of 
the system - wotild have no effect on the flow. Also it could have important consequences 
for related systems such as Rayleigh-Benard convection with through-flow^^’^^. Based on 
numerical solution of the CGL equation it was argued in ref. 15 that thermal noise may 
play an important role in the Taylor-Couette experiments. However, as stressed by the 
authors, their result provides only an order of magnitude estimate since the noise term 
used in the CGL simulation was not derived rigorously from the Navier-Stokes equations. 
We note that thermal noise is also believed to be important in recent experiments in binary 
fluid convection.^® 

In this Letter we will first study, without reference to any particular physical system, 
the CGL equation with a noise term that is delta-correlated in space and time. We will 
present for the first time analytic results for this system under conditions when the system 
is convectively unstable. Since the CGL equation is a generic equation and considering the 
facts that, as noted previously, noise is an element common to all physical systems and all 
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systems with nonzero group velocity are convectively unstable at onset of the instability, 
these results will be applicable to a wide variety of physical systems. 

Next we will focus attention on the particular system mentioned previously (i.e. 
Taylor-Couette flow with am axial through-flow) and rigorously derive the noise term for 
the CGL equation, as well as the coefficients of the CGL equation (including the nonlinear 
term). Based on these restdts, it appears that the noise-sustained structure seen in the 
experiments, although extremely small - being roughly am order of magnitude larger than 
thermal noise - is not of thermal origin. However, even if the structure in these experi- 
ments is not thermal in origin, it seems likely that experiments cam be designed in which 
structure is sustained by thermal noise (i.e. molecular motion), considering the fact that 
the noise in these experiments is already so extremely small. 

The complex Ginzburg-Landau equation^*“^® with a noise term is 

At = taA-vA^-\-hAxx—c\A^A + ^ ( 1 ) 

where o, 6, and c are in general complex coefficients, v is the group velocity, e = R — Rc 
measures the “distance” above onset of the instability (where R and Rc aure the control 
pairameter amd critical value of that parameter, respectively), ^(«,t) is a complex thermal 
noise term, and A{x,t) is the slowly varying complex aunplitude of a plane wave solution 
at criticality. The conditions satisfied by ^ aure (^(a:,<)^*(x',f')) = a^8{x — x')S{t — <'), 
(^(®,<)|(®',t')) = 0, and (^(r,t)) = 0, where * refers to the complex conjugate, 8 is the 
Dirac delta function, and the amplitude of the noise level is ^ven by cr. These condi- 
tions will be satisfied if (^r(a:,f)^r(a:',^')) = (6(®j0^t(a:',f')) = {a^ I2)8[x - x')8{t - <'), 
(^r(®jf)^i(!c',t')) = 0, and (^,.(r,t)) = (^i(r,<)) = 0, where the subscripts r and i refer to 
the real amd imaginairy parts, respectively. As shown by Graham,^® the fluid equations in 
the presence of thermal noise and near onset of the instability may be reduced to eq. (1). 

Assuming that |A| is sufficiently small so that the nonlinear term may be neglected. 
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the covariance function K{xyx\t) ^ (where the angled brackets signify 

the expectation value) satisfies the following equation: 

K-i = 2c(i-pK — v(^Kx “1" -^x') H" “I" b* Kx*x* <7^6(x x ) (2) 


This equation was derived by discretizing the linearized eq. (1) in space, applying the 
results for a set of coupled ordinary differential equation, and then returning to the 
continuous spatial limit . 

Assuming that the system is convectively unstable, i.e. 0 < ea, < v 26 ,/( 4 | 6 | 2 ), 2-3 the 
stationary solution of eq. ( 2 ) for the semi-infinite interval [ 0 , oo) with boundary conditions 
K{x,0,t) = K{0,x',t) = 0 is 


K{x,x') = ® 


7 T-‘ 


' flh fll-' ^ ^ ^ 

Jo Jo + {k — fc')2 7 ^ -f (A: -i- fc')^ / 


(3) 


where 7 = vbr/\b\'^ and a = ^/v'^br/{2\b\^) - 2ear- After taking x = x', making a change 
of variables, and performing one of the integrations, eq. (3) may be reduced to 

s(sx) — ®^cosh(|^sx) 




i : 


cos( 


(4) 


y / s ^ + T }^ (s2 -f 72) 

where rj = (-v/56^ This equation gives the variance K{x,x) = (|A(x,<)p) at the 
point X. For large x eq. (4) may be written as the asymptotic series 

1 


K(x,x) 


tr^j 




y/ 2 TTT} |f>| -rj^ 72-1- 
where the /?^ are the coefficients in the Taylor expansion 
1/2 1 


where 


ho 


(5) 


1/277-1 - 5 V 72 - iq + rjy 72-|- s +(5)2 72-1-3 


hs :)- p -- 
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^±( 9 ) = + 9 ) ± ^\/9(2t7 + 9 ) 
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where r(|) = y/n and F(^ + |) = (£ — ^)F(£ — |), F bding the gamma function. The 
coefficient of x in the exponential, 7 — 7 ;, is precisely twice the spatial growth rate of 
the amplitude of the most rapidly growing mode^’^ as one would expect. Although this 
expression for the variance may look somewhat formidable, note that the leading term in 
the asymptotic series is given simply by replacing the sum in eq.(5) by 1 and, as we shall 
see, this may be all that is necessary in many cases. Another feature that should be noted 
is that the variance does not increase purely exponentially in space for large x , but rather 
as exp ((7 — T})x)ly/x. 

We now focus attention on Taylor-Couette flow with an axial through-flow. This sys- 
tem consists of two concentric cylinders with inner and outer radii rj cmd r„, respectively, 
with the inner cylinder rotating with velocity Wj, and with an imposed through-flow in 
the axial direction. The radial, azimuthal, and axial coordinates are denoted by [r,9,x), 
respectively. In the parameter regime of interest, the flow will be axisymmetric and all 
derivatives with respect to 0 will vanish. To derive the CGL equation the Navier-Stokes 
and continuity equations for the deviation of the velocity and pressure about the station- 
ary background flow are written in cylindrical coordinates and the velocity, distance, time, 
and pressure are scaled with v/d^ d, d^ fv, and Jd?, respectively, where d is the gap dis- 

tance between the cylinders and v is the kinematic viscosity. There are three independent 
parameters for this system: 1) the scaled inner cylinder velocity or azimuthal Reynolds 
number Rcaz = Vid/v, 2) the scaled average axial velocity of the stationary background 
flow or axial Reynolds Rcax = Wdjv, and 3) the scaled inner cylinder radius rj/d, which 
may be written in terms of the radius ratio Tojri as Ti/d = l/[(’’o/’'i) — !]• 

As Rcaz is gradually increased for a given Rcax and r^/d, at Reaz,c the critical value 
of Rtaz, the stationary background flow becomes unstable. We define e that appears in 
the CGL equation (1) as e = Rcaz — R^az,c This measures the “distance” above onset of 
the instability and is used as the expansion parameter in deriving the CGL equation. The 
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linear coefficients a, v, and h are particularly easy to derive, being determined from the 
linear stability problem. Assuming a solution of the linearized fluid equations of the form 
exp{X{k, Reaz)t + ikx), where t and x are the scaled time and axial position, respectively, 
these coefficients are given by a = dX/dRcaz, ^ = —dXi/dk, and 6 = —{lj2)d^Xfdk'^, 
where the derivatives are evaluated at criticality (i.e. at kc and Rea.z,c)- The critical 
vcdues of k and Rcaz are determined from dXr/dk = 0 and Ar = 0. The calculation of the 
coefficient of the nonlinear term, c, is much more involved and we only give a numerical 
value in this paper. For details of the derivation of coefficients of the CGL equation, the 
reader is referred to, for example, ref. 18-20. In terms of the amplitude A appearing in the 
CGL equation (1), the scaled radial, azimuthal, and axial fluid velocities and pressure are 
given by (u,v,uj,p) = A{x,t){U,V,W, P)exp{ikcX — iwct) + c.c., where U{r), V'(r), W{r), 
and P{r) are the radial eigenfunctions of the linear stability problem at criticality, kc and 
Wc are the critical wavenumber and frequency, respectively, and c.c. stands for the complex 
conjugate. Chebychev polynomials were used in the radial direction giving highly accurate 
results. 

In the presence of thermal noise, noise terms (z = r,B,x) must be added to the 
Navier-Stokes equations for the radial, azimuthal, and axial velocity components u, v, and 
w, respectively. From Landau and Lifshitz^^ we find that for thermal noise (writing the 
noise term in cylindrical coordinates) 

A-(0 = 1 ). -f -h {i = r, 0, x) (6) 

where 0), and where the correlation between the random components 

of the stress tensor (in scaled units) are 

I, f)5<‘™>(r', = 
i/^pd [ \ 3/ J ^ 
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where is the Kronecker delta function (i,j = r,^, x), k is Boltzmann’s constant, T 
is the absolute temperature of the fluid, p is the fluid density, and ^ is the ratio of the 
second (or bulk) viscosity to the usual viscosity. Expanding the velocities and random 
components of the stress tensor in e*"*®, where m corresponds to the azimuthal mode 
number, and noting that the m=0 mode dominates, we find — = 1 /( 27 t), in addition 

to all the derivatives with respect to 6 vanishing. 

To find the correlation of the noise term we followed Graham.^® The correlation of 
the noise term is found to be 




27r|Cp 


(g + CKl (3 ” 0^2 + (C - g)^3 + h 


where 


= /* dr- ( -U^- 
J r y r 


dU^ 


dr 


+ 


-J7t 


= I drhl \W^i 


= f dT-( -U^-^+ik^W' + 

J r \ r dr 

U = jdr\ + 


-U'-ik^W' 


r dr 


+ 


-W^ - ^ +ikcU^ 
r dr 


C = J dr {u^u + v^v + w^w) 


( 8 ) 


where Q = kT / (pdi/^) is a dimensionless number related to the thermal noise level and 
where Ut(r), and VTt(r) are the radial agenfunctions of the adjoint stability prob- 

lem at criticality. 

Numerically evaluating the noise term for a radius ratio of ri/r„ = .7376 which cor- 
responds to the experiments of Refs. 13 and 15, <md for an axial Reynolds number of 
Rcax = 3, gives = .32023 Q^(x — x')S{t — t'). For these experiments we also 

have d — .6769 cm, p = 1.04 g/cm®, u = .0158cm^/s, and T = 293 °K which gives 
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Q = 2.30 X 10”^°. The eigenfunctions were normalized such that lJ7(i’/i)| = 1) where 
Tfc = (ri + ro)/2. Specifying how the eigenfunctions are normalized is essential when evalu- 
ating the noise correlation and the nonhnear coefficient, since the values of these quantities 
depend on the normalization of the eigenfunctions. The numerical coefficient of the noise 
correlation is rather insensitive to the value of Rcax small heing, for example, 

.31719 for Rcax = 0. We note that the noise correlation is independent of C, which is prob- 
ably a reflection of the incompressibility of the flow. This is fortunate since the second 
viscosity is unknown for liquids at low frequencies.^* 

Evaluating the linear coefficients for Reax — 3 gives a = .31448 -(- 6.806 x 10 
V — 3.69024, and b = 1.91052 -|- .14070i. The nonlinear coefficient is found to be c = 
.45127 + 4.444 x 10“*i. We also find that fee = 3.13713, Reaz,c = 84.3035, and = 
11.0294. The real parts of the coefficients a, b, and c, and kc and Reaz,a are found to be 
rather insensitive to the value of Rcax for small Rcaxt being, for example, = .31372, 
br = 1.91094, Cr = .45192, kc = 3.13620, and Reaz,c = 84.0148 for Rcax = 0. The 
group velocity, the imaginary parts of the coefficients, and the critical frequency are 0 for 
Rcax — 0 and are approximately proportional to Rcax for small Rtax- 

Figure 1 (solid line) shows the root mean square (RMS) average of the radial velocity 
at r = r;, as a function of x from a numerical simulation of the CGL equation (1) using 
the above values for the coefficients and thermal noise term and using c = 4.668. The 
RMS average of the radial velocity at r = rji is related to the RMS average of \A\ by 
(u(r/i,x,f)^)*/^ = V^(|j4(x,t)p)^/^ by virtue of the normalization taken for the radial 
eigenfunctions (see discussion following eq. (8)). As can be seen, the noise is spatially 
amplified resulting in a noise-sustained structure. In the absence of a continuous source of 
noise, the solution is zero everywhere. Figure 1 (dashed line) shows at r = as 

given by the leading term in eq. (5). As cau be seen the agreement is excellent between 
the analytic result and the numerical solution of the CGL equation for small |A|. In order 
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to see the range of x for which eq. (5) is valid, fig. 2 shows from the numerical 

simulation (solid line), from the leading term in eq. (5) (dashed line), and from eq. (5) 
keeping three terms in the asymptotic series (ie. n = 2) (dotted line). Although the RMS 
average as given by the leading term in eq. (5) deviates somewhat from the numerical 
solution for smaller x, the RMS average as given by the numerical solution and the RMS 
average as given by eq. (5) with n = 2 show excellent agreement along most of the curve. 
Referring back to fig. 1 we see that the value of x at half maximum is about 26.5. For 
the experiments of Refs. 13 and 15, this distance is about 21.^^ Referring to eq. (5) and 
noting that {-y — r))/2 = .5625, this implies that the noise amplitude in the experiments is 
about (21/26.5)^ /^exp( .5625 x (26.5 — 21)) rs 20 times thermal noise. 

In conclusion, we have studied the complex Ginzburg-Landau equation with a thermal 
noise term under conditions when the system is convectively unstable. Analytic results were 
derived for the variance. Since the CGL equation is a generic equation, and considering 
the facts that noise is an element common to all physical systems and all systems with 
nonzero group velocity axe convectively unstable at onset of the instability, these results 
will be applicable to a wide variety of physical systems. The coefficients and thermal 
noise term for the CGL equation were determined for Taylor-Couette flow with an axial 
through-flow and comparison was made to experiment. Although the effective noise level 
in the experiments of Ref. 13 and 15 appears not to be thermal in origin, the noise level 
is nonetheless extremely small, being roughly an order of magnitude larger than thermal 
noise. Therefore, even if the structures in the experiments are not thermally sustained, it 
seems likely that an experiment can be designed in which thermally-sustained structure 
does exist. 

I thank Wai-Ming To for providing a code for doing the linear stability analysis and 
for his reading of and helpful comments on the paper. I thank Ken Babcock for stimulating 
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Fig. 1. — RMS average of the radial velocity at a radius midway between the cylinders from a numerical simulation ( 
CGL equation (solid line) and from the leading term (i.e., n = 0) in the asymptotic series Eq. (5) (dashed line). 
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